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1. Introduction 


With the acceleration of urbanization and improvement of 
people’ living standard, a larger proportion of building energy 
— _ , consumption will be needed to keep a comfortable indoor envir- 
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Nomenclature 

A specific surface area per unit volume [m7 !] 
Cp specific heat [J kg~! K7!] 

Coat saturation specific heat [J kg~' K~1] 

dh hydraulic diameter [m] 

Dis mass diffusion coefficient [m? s7 !] 

Di thermal diffusivity [m? s71] 

g gravity [ms~7] 

G specific mass flow rate [kg m~? s7 !] 

G’ mass flow rate [kg s7 !] 

h specific enthalpy [kJ/kg] 

he enthalpy of humid air in equilibrium with liquid 


desiccant [kJ kg] 
Neeff effective saturation enthalpy [kJ kg], Eq. (18) 
height of the dehumidifier [m] 


k heat conduction coefficient [W m7! K~'] 

L length of the dehumidifier [m] 

Le Lewis number 

m water condensation rate [gs7!] or per unit cross- 
sectional area [g m7? s~ 1], Eq. (23) 

m* capacity ratio similar to the one used in sensible heat 
exchangers 

M molecular weight [g mole~ !] 

NTU number of transfer units 

Nu Nusselt number (dimensionless) 

P pressure [Pa] 

Pi partial vapor pressure in air [Pa] 

P; partial vapor pressure over the solution [Pa] 

P: total pressure [Pa] 

P dimensionless vapor pressure difference ratio 

Q heat transferred from solution to water [kW m7 !] 

T temperature [K] 

T dimensionless temperature difference ratio 

u velocity [m s71] 

V volume [m?] 

w humidity ratio [kg H20 kg~! dry air] or the width the 


dehumidifier [m] 


surface for breeding mildew and bacteria and so on [2]. Thus, to 
reduce the energy consumption in buildings and improve the 
indoor air quality, the liquid desiccant assisted air conditioning 
system has drawn more and more attention [3-7]. 

The major component of interest regarding heat and mass 
transfer of such a system is the dehumidifier. Compared with the 
experimental research, the simulation method is more time and 
cost saving. Also, some parameters in the dehumidifier interior can 
be observed by simulation while it is impossible to be achieved by 
experiment. Most importantly, the verified simulation models are 
effective tools to assess and optimize similar dehumidifiers. There- 
fore, a large amount of studies have been done to establish 
reasonable mathematical models for evaluating the liquid desic- 
cant dehumidifiers. However, there is short of detailed and specific 
summary of the models until now. Thus, it is meaningful to classify 
and assess the models, which will provide useful suggestion for 
future research. 

In the paper, the function principle of the liquid desiccant 
dehumidifier is introduced firstly. Based on whether there is heat 
removal, the dehumidifier is divided into two types: adiabatic and 
internally cooled dehumidifier. Correspondingly, the models are 
summarized in two respects in terms of the different structures. 
For each model, the assumptions, governing equations, boundary 


We humidity ratio of humid air in equilibrium with liquid 
desiccant [kmol H20 (kmol air)~'] 

We eff effective humidity ratio [kmol H20 (kmol air)~*], 
Eq: (19) 

X desiccant concentration [kg desiccant kg~ solution] 

Xw concentration of water in solution 
[kg water kg~'solution] 

Xy concentration of water vapor in air [kgH20 kg- !] 

Greek letters 

ac heat transfer coefficient [W m7? K7 !] 

aç heat transfer coefficient corrected for simultaneous 
mass and heat transfer [W m7? K7!] 

ap mass transfer coefficient [kg m7? S71] 

ap molar mass transfer coefficient [kmol] m7? S71] 

9 dimensionless temperatures (T—T,)/h, Eq. (33) 

p density [kg m~?] 

u dynamic viscosity [kg m~!s~1] 

v kinetic viscosity [m° s71] 

E air side effectiveness 

EHE heat exchanger effectiveness 

À latent heat of vaporization [kJ kg~!] 

ô film thickness [m] 

A change of or difference between parameters 

Subscripts 

a air 

c critical 

f cooling fluid, like water, air, refrigerant 

i inlet 

int interface 

o outlet 

p primary air 

r secondary (return air) 

s desiccant solution 

v water vapor 


conditions and other relevant information are provided. The 
applied conditions, development process, and research status of 
the simulation models are also presented. In addition, some 
suggestions are put forward for the model improvement. 


2. Problem description 


2.1. The mechanism of heat and mass transfer in liquid desiccant 
dehumidifier 


It is well know that in the dehumidifier, complicated heat and 
mass transfer occurs. The driving force for heat transfer is the 
temperature difference between the air and desiccant solution, 
and for mass transfer is the water vapor pressure difference 
between the air and the surface of the desiccant solution. The 
most classic and typical mass transfer theories include film theory, 
penetration theory and surface renewal theory. The theory used 
most for the dehumidifier is the film theory. It is Nernset [8] who 
proposed the film theory first in 1904. He assumed that the whole 
resistance of mass transfer in a given phase lied in a thin and 
stagnant region of that phase at the interface. This region is called 
film. Based on it, Whiteman [9] developed the two-film theory. 
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The specific transport mechanism is shown in Fig. 1. Where, Pg is 
the partial pressure of component B in the gas phase and Xz is the 
mole fraction of component B in the liquid phase. 

The two film theory is very easy to understand and apply, but it 
has several drawbacks. Firstly, it is not reasonable to predict that 
the rate of mass transfer is directly proportional to the molecular 
diffusivity. Secondly, it is difficult to decide the thickness of the 
two laminar sub layers by experiment. Finally, the convective mass 
transfer in the thin films is neglected, so the theory is only suitable 
for the steady mass transfer process. 

In the dehumidification process, some quantity of heat is given 
out as well, including the phase change heat and dilution heat. In 
all of the dehumidifier models, the dilution heat is neglected as it 
is much smaller compared with the phase change heat of water 
vapor [10]. 


2.2. Structure of the dehumidifier 


As for the structure of the dehumidifier, according to whether 
there is heat output, the dehumidifiers can be classified into 
adiabatic and internally cooled dehumidifier [11]. The diagrams 
of two dehumidifier structures are shown in Fig. 2. 

In an adiabatic dehumidifier, the air and liquid desiccant 
contact directly with each other. In the early stage, the research 
was concentrated on the structure of spray tower as a result of its 
simple construction and large specific surface area [12]. However, 
in the spray tower, the desiccant solution is generally broken into 
small droplets, so the problem is sometimes serious of mist 
generation and carryover of liquid droplets in the air stream. Then, 
the packed tower was used widely because it is more compact and 
can provide a higher residence time, lower liquid pressure loss and 


Interface 


Mass Transfer Direction 


Fig. 1. Schematic diagram of two-film theory. 


Concentrated solution 


Packing 


Diluted solution 


lessen the carryover problem. In 1980, Factor and Grossman 
verified the possibility of employing the packed tower as dehu- 
midifier by theoretical analysis and experiment [13]. As for the 
padding materials, the random packing like ceramic [|14], plastic 
and polypropylene pall ring are popular first [15,16]. Then some 
structured packing materials are employed to optimize the flow 
and reduce the resistance in the dehumidifier, such as the stainless 
steel corrugated orifice plate [17], celdek [18,19] and so on. 

In the adiabatic dehumidifier, the temperature rise of the 
desiccant, resulted from the latent heat, worsens its dehumidifica- 
tion performance, thus its dehumidification efficiency is relatively 
lower. A solution is to increase the desiccant flow rates to achieve 
good dehumidification levels. However, the high desiccant flow 
rates and the followed higher flow rates of the regenerated 
desiccant solution reduce the coefficient of performance of the 
liquid desiccant cycle [20,21]. In addition, the desiccant particles 
are much easier to be entrained by the air and therefore pollute 
the indoor environment. 

To solve the above problems, the internally cooled dehumidifier 
was developed. In an internally cooled dehumidifier, besides the 
contact between air and desiccant, some cold source which can 
provide cool fluid like air or water is added to take away the latent 
heat produced in the process of dehumidification, which can be 
regarded as an isothermal process in general. The internally cooled 
dehumidifier has been popular since the 1990s [22,23]. As the 
latent heat is removed from the dehumidifier, it reduces the 
temperature rise of the solution and air, resulting in efficiency 
improvement [24]. Meanwhile, it allows lower desiccant flow rates 
in the internally cooled dehumidifier so as to reduce the pollution 
problem. But the internally cooled dehumidifier has more com- 
plicated structure than the adiabatic one. For example, to increase 
the contact area, a fin structure is widely used in the internally 
cooled dehumidifier or other heat and mass transfer devices 
[25-29]. 


3. Models for adiabatic dehumidifier 


There are mainly three types of mathematical models, includ- 
ing the finite difference model, effectiveness NTU (e—NTU) model 
and the simplified solutions. 


3.1. Finite difference model 


In 1980, Factor et al. [13] promoted a theoretical model to 
predict the performance of a countercurrent packed column air- 
liquid contractor, based on the model for adiabatic gas absorption 
put forward by Treybal in 1969. In the model, the whole dehumi- 
difier is divided into n parts, as shown in Fig. 3. 

To simplify the complexity of the heat and mass transfer 
process, several assumptions were made: (1) the flow of air and 
desiccant were assumed as the slug flow, (2) the process was 
adiabatic, (3) the properties of the gas and liquid were assumed 


b 


Dry Air 
fom Concentrated solution 
ye 
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Packing 
fluid 
Moist Air Lae Dilated solution 


Fig. 2. The structure diagram of two dehumidifiers. 
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constant across the differential element, which meant the gradi- 
ents only exist at the z direction, (4) both of the heat and moisture 
transfer areas were equal to the specific surface area of the 
packing, (5) it was negligible of the non-uniformity of the air 
and solution flows, (6) in the flow direction, no heat and moisture 


Fig. 3. Heat and moisture exchange model in the countercurrent adiabatic 
dehumidifier. 


transfer occurred, (7) the resistance to heat transfer in the liquid 
phase was negligible, and (8) the interface temperature was equal 
to the bulk liquid temperature. Based on the above assumptions, 
the main governing equations were stated. 

According to the mass balance in the control volume, 


dG; = GadW (1) 


According to the interface mass and sensible heat transfer rates, 
the air humidity change was, 


dW _apMvA n (: AP: R a 


dz Ga 1— Pa/Pt 


According to the interface sensible heat transfer from the air to 
solution side and the energy balance on the gas side, the air 
temperature change was, 


dTa a gA(Ta —Ts) 
dz GaCpa (3) 
aA — GaCpv(dW /dz) m 


€ T= eXp[GaCp y(dW /dz)/(acaA)] 


where acq and a, are the heat transfer coefficient (sensible) of 
the gas side and that coefficient corrected for simultaneous mass 
and heat transfer by applying the Ackermann correction, which is 
one method to take into account the effect of mass transfer on the 
temperature profile with an Ackermann correction factor. 

Finally, the boundary conditions were: z=0 T;=Tsj, Gs= Gsi, 
X Xi 2 H. Ta Tai Ga Gai, W W;. 

Since the above differential equations cannot be solved analy- 
tically, the most basic solution is a numerical integration along the 
height of the dehumidifier. To begin the calculation, one end of the 
dehumidifier must be chosen as the start point. For the counter- 
current flow pattern, it needs to presume the outlet conditions for 
one of the fluids. Solving the above equations from the bottom to 
the top of the dehumidifier, with the boundary conditions, a set of 
calculated inlet solution parameters are obtained. By comparing 
the calculation results with the real values, the supposed existing 


Input the inlet 
parameters of the 
air and solution 


$ 


Suppose the outlet parameters of 
the solution and initialize the 
whole flow field 


l 


Input the physical parameters of > 
the air and solution 


With the governing equations, 
calculate the relevant parameters 
from j=1 to j=n 


Iteration 


Obtain the outlet parameters of 
the solution 


Do not meet the 


Compare the calculation values of the 
outlet solution with the real values 


accuracy requirements 


Meet the accuracy 
requirements 


The supposed values of the 
outlet parameters are correct 


Į 


Calculate the outlet parameters 
of the air 


Fig. 4. Calculation flow chart of countercurrent pattern. 


Y. Luo et al. / Renewable and Sustainable Energy Reviews 31 (2014) 587-599 591 


L zZ 
a Desiccant 
L 


Air 


X 


Fig. 6. A two dimensional schematic of the cross flow dehumidifier/regenerator 
[39]. 


solution variables are adjusted. The calculation will last until the 
final results are very close to the real values. And the general 
calculation flowchart for the finite difference model of the counter 
flow pattern can be summarized in Fig. 4. 

In the later study, Gandhidasan et al. [30] utilized the similar 
model to study various parameters on the packing height of 
the packed tower. Then, Lazzarin et al. [31] gave more specific 
explanation of the calculation method in Appendix A of the litera- 
ture. Oberg and Goswami [32] applied a finite difference model 
similar to Factor and Grossman's to verify the experimental results. 
Taking account the insufficiently wetted packing and the different 
factor when transfer the k-type mass transfer coefficients to the F- 
type one, Fumo and Goswami [33] improved Oberg and Goswami’s 
mathematical model by modifying the transfer surface. In addition, a 
correction factor CF was introduced to modify the correlation of the 
wetting surface. By comparing the results of simulation and experi- 
ment, it was found that the calculation results of the adapted model 
agreed well with experimental results. 

All of the above models introduced a coefficient aç to describe 
the simultaneous heat and mass transfer with the Ackermann 
correction. Khan and Ball [34] promoted another solution to deal 
with the simultaneous transfer process. Both heat and mass 
transfer processes were assumed to be gas controlled, so the 
interface temperature was the temperature of the bulk liquid 
and the heat transfer rate across the air film from the bulk air to 
the interface was equal to that entering the liquid side, 


GaCp,adTa = ac aA(Ts —Ta)dz (5) 


Similarly, the mass transfer across the interface was equal to 
the change in humidity ratio, 


GadW = apaA(We — W)dz (6) 
Then, the humid air specific enthalpy change could be written 

as, 

dha = Cp,adTa +dW ` [Cp (Ta = T;)+4] (7) 


By substituting Eqs. (6) and (7) to (5), the air enthalpy change 
along the flow direction was obtained. Here it is rewritten in a 
simpler way appeared in the later literature, 


ðhħa NTU-Le 


1 
re N E (e-noa) we- w] ® 


In the above equations, Le and NTU were defined as 


ac 


Le= 9 
AOR (9) 
NTU = (10) 


a 


In this way, the coupled heat and mass transfer were consid- 
ered together. In the following research, this handling method is 
more popular than the Ackermann correction. 

The finite difference model has been widely used for the 
countercurrent dehumidifier [35-38]. For cross flow configuration, 
Liu et al. [39] proposed a model for the heat and mass transfer 
process in a cross flow adiabatic liquid desiccant dehumidifier/ 
regenerator. The physical and mathematical models are described 
in Figs. 5 and 6, respectively. 

The governing equations of energy, water content, and solute 
mass balances in a differential element were, 


Gi. dha 1 Gs) _ 


H oz L Ox E aD 
G, oW 1 0G, _ 

H z L” a2) 
d(G; :-X)=0 (13) 


The energy and mass transfer in the interface of the air and 
desiccant solution were expressed in Eq. (8) and the following 
Eq. (14), 

W NTU. (We-W) (14) 

Like some other papers, Le was supposed to be one in the 
model. However, the value of NTU was correlated based on the 
corresponding experimental data in the paper. 

Niu [40] also established a two dimensional mathematical 
model for the cross-flow adiabatic dehumidifier, and the mass 
transfer coefficient was gained from the experimental data. Woods 
and Kozubal [41] applied the finite difference model to study the 
performance of a desiccant-enhanced evaporative air conditioner 
and the simulation results showed good agreement with the 
experimental ones. 


3.2. Effectiveness NTU (e—NTU) model 


Stevens et al. [42] reported an effective model for liquid- 
desiccant heat and mass exchanger, which was developed from a 
computationally simple effectiveness model for cooling towers 
[43]. Except for the assumptions of the finite difference model, two 
additional assumptions were included. One was the assumption of 
the linear relationship of saturation enthalpy and temperature, the 
other one was the neglect of the water loss term for the solution 
energy balance. In addition, an ‘effective’ heat and mass transfer 
process was assumed. 
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The main equations and calculation process of e-NTU model 
were summarized as follows, 


(1) Calculated the Number of Transfer Units (NTU) by Eq. (10). 

(2) In terms of the similarity with the heat exchanger, the 
effectiveness of the countercurrent flow dehumidifier could 
be expressed as 


1 —e-NTUA —m*) 


= Time NUO m (15) 
where m* was a capacitance ratio, defined analogous to the 
capacitance ratio used in sensible heat exchangers, and it was 


given in the following equations, 


GC. 

m* = =3 (16) 
G; iCp,s 

where Csa was the saturation specific heat, and 


Csat = (dhe /dT;). 
(3) With NTU and e, the air outlet enthalpy could be obtained with 
the following equation, 


Nao = hai +e(lhe— Nai) (17) 


(4) Used an energy balance to calculate the solution outlet 
enthalpy. 
(5) Then an ‘effective’ saturation enthalpy was found by 
Nao — Nai 


heeff = hai t7 S=nTU (18) 


(6) Using the enthalpy and saturated conditions, the effective 
humidity ratio Ye er; could be obtained. 

(7) Then, by the following equation, the air outlet humidity ratio 
could be calculated, 


Wo = We eff +(Wi-Weerpe “NTT (19) 


(8) With the mass balance and known inlet and outlet parameters, 
all of the outlet parameters were acquired. 


In the later study, Sadasivam and Balakrishnan [44] pointed out 
that the definition of NTU based on the gas mass velocity in 
Stevens’s model was not appropriate when the minimum flow 
capacity was the liquid [45]. Thus, the gas mass velocity G, in 
Eq. (10) was changed to the minimal mass velocity of gas and 
liquid. 

As for the e-NTU model, there are much fewer literatures than 
the finite difference model. In the following study, Ren [46] 
developed the analytical expressions for the e-NTU model with 
perturbation technique. The model accounted for the nonlinea- 
rities of air humidity ratio and enthalpy in equilibrium with 
solutions, the water loss of evaporation and the variation of the 
solution specific heat capacity. 


3.3. The simplified models 


It can be found that both the finite difference model and e—NTU 
model require numerical and iterative computations. Thus, both of 
them are not suitable for hourly performance evaluation. 

Khan and Ball [34] developed a simplified algebraic model. 
With the finite difference model, about 1700 groups of data were 
analyzed. The following functions were deduced, 


Wo = Not+n W; +n Tsi + n3Te; (20) 


Wo = Mo +m Wi +mM2Tao+M3Te (21) 


The above two equations can be employed to predict the exit 
air temperature and humidity ratio easily. However, as they were 
fitted based on the data with certain operating solution concen- 
tration and solution to air mass ratio, they might not be suitable 
for other conditions. 

Liu et al. [47] fitted some empirical correlations to estimate 
the performance of a cross-flow or counter-flow liquid desiccant 
dehumidifier. The essence of the methodology is obtaining the 
empirical expression of enthalpy and moisture effectiveness by 
experiment. 

Gandhidasan [48] reported a very simple analytical solution to 
predicate the rate of moisture removal for the dehumidification 
process. Referred to previous work [49], the author promoted 
the dimensionless moisture and temperature difference ratios. 
By combining the aforementioned two definitions, the energy 
balance equation could be expressed as follows, 


My A 
Ma Pt 
According to the literature [50], the relationship between the 
rate of moisture removal m and the partial pressure of water vapor 
could be deduced as 

mP: Ma 

Psi = Pai -> r 

S, a,l G-P= M, 
In addition, the desiccant outlet temperature was easily calcu- 

lated, as shown in the equation, 


= = G 
Cpal (Tai Ti) t P(Pai P34) = c CpslTs.o Tsi) (22) 
a 


(23) 


T. = Isi EHET ci 
60 (ene) 


Finally, substituting Eqs. (23) and (24) to (22), the moisture 
removal rate m was given as, 


m= T G;Cp,sEHE 
A| — ene) 


The method was rather simple yet it involved lots of assump- 
tions and limitations. Except some common assumptions, it also 
required that the desiccant inlet temperature was different from 
the air inlet temperature, and the desiccant temperature leaving 
the dehumidifier was equal to that of the desiccant entering the 
heat exchanger. 

Chen et al. [51] constructed an analytical model on the basis of 
the finite difference model for countercurrent and concurrent flow 
pattern. The physical model is similar to that of Fig. 3. Firstly, a 
mathematical model was built following the model promoted by 
Khan and Ball [34]. 

Then, two parameters were introduced for derivation conve- 
nience, 


Ka=le Cpa’ Tata W (26) 


24) 


ee G,CpaT Tas—T30) (25) 


Ke =Le: Cpa ` Ts+4:We (27) 


By combining the mass, energy conservation equations, mass 
and energy transfer equations at the interface, the change of Ke 
along the flow direction was, 
dKe m*NTU 


1 
dz = H haitz Ke Keo) ł (Le 1)Cpa ‘TaK. (28) 


By integrating (28) from 0 to H or z, Ke at the outlet and along 
the z-axis were acquired. Based on the above results, the distribu- 
tion of air enthalpy, air humidity, and temperature in a counter- 
current adiabatic dehumidifier were obtained. Then, the 
distribution of the solution parameters can be calculated. In the 
paper, the analytic solution of the concurrent flow heat and mass 
transfer process were also given out. 

Ren et al. [52] derived a new analytical solution from one- 
dimensional differential model. By introducing some 
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Table 1 
Comparison of the mathematical models for adiabatic dehumidifier. 


Classification Assumption Iteration 
Finite difference model Least Extensive 
e-NTU model More Less 
Simplified models Most No 


Table 2 
Detail information of the mathematical models for adiabatic dehumidifier. 


Classification Model Flow pattern 
Finite difference model Factor et al. [13] Counter 
Oberg et al. [32] Counter 
Fumo et al. [33] Counter 
Khan et al. [34] Counter 
Liu et al. [39] Cross 
e-NTU model Stevens et al. [42] Counter 
Classification Model Flow pattern 
Simplified models Khan et al. [34] Counter 
Liu et al. [47] Counter or cross 
Gandhidasan [48] Counter 


Chen et al. [51] 


Ren et al. [52] Counter 
Babakhani et al. [54] Counter 
Liu et al. [56] Cross 


dimensionless and dimensional groups, the conventional equa- 
tions of one-dimensional model were transferred to two coupled 
ordinary differential equations, whose general solution are as 
follows, 


AWu=Ci e^ NTU: + Co e2NTUz (29) 


Ad = —KyCye® TY: 4 Ky Ca e2NTUe (30) 


The model is just suitable for the case where the solution flow 
rate and concentration change slightly as it assumed that the 
variation of the equilibrium humidity ratio of solution depended 
only on the change of the solution temperature. 

Babakhani and Soleymani developed analytical models for the 
counter flow adiabatic regenerator [53] and dehumidifier [54]. The 
analytical solution was deduced on the basis of the differential 
equations, including the mass balance, air humidity and tempera- 
ture change equations listed in (1)-(3), and liquid desiccant 
temperature and concentration change derived from the mass 
and energy conservation equations. To achieve the simplification, 
two main assumptions were applied, including the assumptions of 
the dilute gas phase and constant equilibrium humidity ratio on 
the interface. Then the above equations could be solved and the 
integrated analytical solution were, 


W = Win +(Wi- Wint)exp(— aMNTU;) (31) 
p — 
Ta = C1 + C2€X ONTU — — ex MNTU 32 
a 1 2exp( z2) (aM) — ao p(— a z) (32) 
1 aMB = 

T= Rale! — C20exp( — ONTU;) + Gi? aMo exp(—aMNTU;)]+Ta 
(33) 
InX = — Rm(Wi— Wint)eXp(— aMNTU,) +C3 (34) 
Gs = Ga(Wj — Wine)exp(— aAMNTU;)+ C4 (35) 


Counter or concurrent 


Accuracy Applied situation 

Best Component design and operation optimization 
Better Component design and operation optimization 
Worst Annual assessment 


Dimensionality Treatment of coupled, heat and mass transfer 
Ackermann correction 

Introduction of enthalpy 

Ackermann correction 

Ackermann correction 

Introduction of enthalpy 


One-dimensional 
One-dimensional 
One-dimensional 
One-dimensional 
Two-dimensional 


One-dimensional Introduction of enthalpy 


Simplified method 


Correlations based on simulation results 

Correlations based on experimental results 

Introduction of dimensionless parameters 

Introduction of parameters similar to air enthalpy 
Introduction of dimensionless parameters 

Transfer the differential equations to nonlinear equations 
Similar method of cross-flow heat exchanger 


With the above solutions, the profiles of the outlet solution and air 
parameters were available. 

Based on the above model, Babakhani [55] also developed 
another analytical model which was well suited to the high 
desiccant flow rate conditions. 

Liu et al. [56] developed an analytical solution for a similar 
cross-flow packed bed liquid desiccant air dehumidifier, whose 
numerical model had been reported in literature [39]. In the 
present work, it was regarded that the desiccant mass flow rate 
and concentration kept constant in the whole process. Thus the 
Eqs. (12)-(14) were got rid of, and only Eqs. (8) and (11) were left 
for calculation, which were rewritten as 


alta Hale _ 
ier ET rami G0) 
aha NTU 

op LT Mahe) (37) 


It was found the above control equations had high similarity 
with those of the cross-flow heat exchanger. Thus the methods in 
literature [57,58] were referred to. As the solution expressions 
were a litter complicated, here they will not be presented. 

Recently, Wang et al. [59] developed a simplified yet accurate 
model for real-time performance optimization. Levenberg—Mar- 
quardt method was used to identify the parameters. The proposed 
model is suggested to be employed in real-time performance 
monitoring, control and optimization. Park and Jeong [60] also 
developed a simplified second-order equation model as a function 
of operation parameters to study their impact on the dehumidi- 
fication effectiveness. 


3.4. Summary 


To sum up, there are three models to simulate the adiabatic 
dehumidifiers. The general comparison of them is presented in 
Table 1. As a matter of convenience, the information of some 
representative models is summarized in detail in Table 2. 
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The finite difference model is used most frequently for its 
accuracy. However, it involves complicated iterative process, so it 
is only suitable for the component design and operation optimiza- 
tion. The common assumptions have been stated in Factor’s model. 
However, some other simplifications or improvements are made to 
satisfy the real conditions. It can also be observed that counter 
flow configuration is the most commonly used flow pattern, and 
it can be described by the one-dimensional model. For the 
dehumidifier with cross flow configuration, the problem is gen- 
erally solved by the two-dimensional model. 

For the e~NTU model, two additional assumptions are included 
as mentioned above. Thus, the model is less accurate than the 
finite difference model. However, compared with the finite differ- 
ence model, the calculations are dramatically less extensive. Thus, 
the e-NTU model has great potential for saving time. But with the 
development of the computer technology, the calculated amount 
of the finite difference model can be accepted due to its accuracy. 
Therefore, much fewer researches have been done on the -NTU 
model in the later studies. 

From Table 2, it is concluded that different dimensionless 
parameters have been introduced in the deviation processes of 
the simplified models. Meanwhile, the additive assumptions are 
needed for simplification. Therefore, they are only applicable for 
certain operation conditions, and different models can be chosen 
to be suited to the real situation. The biggest advantage of these 
models is that the iteration is avoided. As a result of their high 
efficiency, they are often used to predict the annual energy 
consumption of an air-conditioning system. 


4. Models for internally cooled dehumidifier 


In the internally cooled dehumidifier, some cooling fluid is 
introduced to remove the heat produced by the vapor condensa- 
tion, as shown in Fig. 7 [61]. Most of the models for the internally 
cooled dehumidifier are developed upon the finite difference 
methods used for the adiabatic dehumidifier. The difference is 
the additional consider of the heat transfer brought by the cooling 
media. In addition, because of the relatively low solution flow rate 
in the internally cooled dehumidifier, it is easier for the solution to 
form a thin film on the surface of the padding wall or plate. Briefly, 
there are also three types of models for the internally cooled 
dehumidifier. 


Fluid in Solution in Air out 
Gr nv; G, > Ts X Ga > TatdT, 
h,+dh, 


W+dw 


/ 


Fluid out Solutionout Air in 

Ge, TtdT¢ G,+dG, Ga, Ta 
T;+dT, ha, W 
X+dX 


Fig. 7. Heat and moisture exchange model in the internally cooled dehumidifier 
[61]. 


4.1. Models without considering liquid film thickness 


The first type ignores the liquid film thickness. Khan and 
Martinez [62] developed a mathematical model to predict the 
performance of a liquid desiccant absorber integrating indirect 
evaporative cooling to achieve an almost isothermal operation. 
The processed air and the solution flowed in countercurrent 
direction while the solution and water flowed in parallel direction. 


a 
Secondary 
air stream 
Primary air 
stream 
Water spray Liquid desiccant spray 
N A 
i \ 
\ 
Primary Secondary 
air stream air stream 


Fig. 8. Schematic diagram of the cross-flow type plate heat exchanger [61]. 


a 


Centreline 


Desiccant 
solution 


Primary air 


Centreline 


dy Thin Water 


Secondary air 


Fig. 9. Schematic diagram of the control volumes considered: (a) primary air- 
solution; (b) secondary air-water [63]. 
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Saman and Alizadeh [63] also established a similar model for a 
cross-flow type plate heat exchanger (PHE) serving as internally 
cooled dehumidifier. The schematic diagrams of the PHE are 
presented in Figs. 8 and 9. It is obviously that it has the same 
principle with that in literature [62]. The only difference lies in the 
configuration difference. Here, Saman’s model will be explained 
in detail as it is more representative. In the paper, the PHE was 
divided into a certain number of control volumes in two orthogo- 
nal directions. Several assumptions were set, including no heat 
transfer with the environment, negligible temperature gradient 
between the solution film and water film, and fully cover of 
solution and water on the plate. 

Thus, the change of the enthalpy and humidity ratio of the 
primary air were, 


dh NTU; - Le 1 

LA i [hes hap)+ a(i 1) ‘Wes-Wp)] (38) 
dW, NTU 

TA =-H s (Wes Wp) (39) 


Also, the mass and energy conservation equations in the 
control volume were given as 


, dW, dG, = 

Cap dx dx (40) 
, dh d(G.h 

Gap _ aa J4Q=0 (41) 


Similarly, the change of the enthalpy and humidity ratio of the 
secondary air, and mass and energy conservation equations were 
expressed as follows 


dh NTU; - Le 1 

ar NTU Le hg, hap a( pe 1) -(Wey—Wo] (42) 
dW, NTU 

vi g Wes W,) (43) 
, dw dG; 

"w dy = ve 


ae dhar d(G;hp) 
ar dy dy 


On the basis of the above governing equations, the discretization 
equations were derived for each control volume. The calculation 
method was similar to the flow chart of the finite difference model 
presented in Fig. 4, while it was more complicated as a result of 
another iteration started by the advance assumption of the cooling 
water temperature. In the paper, it is noted that the film thickness was 
mentioned, but the simulation did not take it into consideration. 

Liu et al. [21] compared the performances of internally cooled 
dehumidifiers with various flow patterns. A representative heat 
and mass transfer model was selected for detail description, as 
shown in Fig. 10. Referred to literature [34], the heat and humidity 
changes of the air were almost the same with the Eqs. (6) and (8). 

The energy conservation equation was expressed as 


g Aha _ 2(G;hs) LOT 
a a Ox Hoy 


Q=0 (45) 


F CpfGFf (46) 

As the mass conservation equation has been explained before, 
here it will not be presented for limited space. And the heat 
transfer between the desiccant solution and the cooling water was, 
oTs NTU; 
y L (Ts—Ty) (47) 

Combined the above equations with the inlet conditions, the 
distribution of the parameters of the air, desiccant solution, and 
the cooling water could be obtained. The model was applied for 


Desiccant in Air out 
(T,, X) SW) 
-m m = - -> 
Water in Water out 
(T) (TdT) 
---- --> 
Desiccant out Air in 
(T,+dT, , X+dX) (T,+dT, , W+dW) 
mm) Liquid desiccant 


— > Humidair 
===> Cooling water 


Fig. 10. Schematic diagram of the control volume for a counter-flow internally 
cooled dehumidifier [21]. 


analyzing the role of flow patterns in depth, which was seldom 
reported in previous literatures. In the later study, Zhang et al. [64] 
employed the above model to predict the performance of an 
internally-cooled dehumidifier. 

Yin et al. [25] built a uniform mathematical model for an 
internally cooled/heated dehumidifier/regenerator which was 
made up of a plate heat exchanger. It is important to point out 
that the author took the non-wetted area into consideration 
by introducing the wetness coefficient. Meanwhile, in a control 
volume, the transfer process in the channel width was considered 
to be symmetrical. In addition, to improve the accuracy of the 
model, the author also applied the correlation of the mass transfer 
coefficient fitted out by experiment [24]. 

Based on the one-dimensional differential equations for the 
heat and mass transfer processes with parallel or counterflow 
configurations, Ren et al. [61] developed an analytical model for 
internally cooled or heated liquid desiccant-air contact units. To 
increase the accuracy, the model took the effects of solution film 
heat and mass transfer resistances, the variations of solution mass 
flow rate, non-unity values of Lewis factor and incomplete surface 
wetting conditions into consideration. 

Recently, Qi et al. [65,66] developed a simplified model to 
predict the performance of the internally cooled/heated liquid 
desiccant dehumidification system. Compared to previous study, 
all of the outlet parameters could be obtained by the correlations 
quickly and accurately. Thus, it is very useful for researching the 
dynamic operation performance of the internally cooled/heated 
liquid desiccant dehumidification system. 

Khan and Sulsona [28] developed a two-dimensional and steady- 
state model for a vapour compression/liquid desiccant hybrid cooling 
and dehumidification absorber made up of the tube-fin exchanger. To 
simplify the model, several assumptions were made reasonably. The 
most critical one was to consider the air and refrigerant flow in 
countercurrent due to the large size in the air flow direction. In this 
way, the complicated problem was simplified to a two-dimensional 
one. The governing equations are very similar to those of literature 
[21] except for the mass conservation equation. There were three 
iterations: one for the real refrigerant exit enthalpy, one for the correct 
local coil surface temperature and the final one for the actual solution 
temperature. 


4.2. Models considering uniform liquid film thickness 


In the second model, the film is considered as a uniformly 
distributed desiccant film. 
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Park et al. [26] developed a three dimensional numerical model 
for simulating the coupled heat and mass transfer in a cross-flow 
internally cooled/heated dehumidifier/regenerator. The schematic 
diagram was presented in Fig. 11. Some assumptions were used 
before listing out the governing equations: the flow was consid- 
ered as laminar and steady; the physical properties for both 
solution and air were constant; species thermo-diffusion and 
diffusion-thermo effects were negligible and thermodynamic 
equilibrium existed at the solution-air interface. Because of the 
relatively small absorption, the TEG solution film thickness was 
simplified as constant and the film mean velocity was also 
unchanged. Also, the velocity gradient in the liquid solution film 
at the interface was regarded as zero and the flow of the liquid 
solution and air was supposed to be fully developed at the start 
place, as presented in Fig. 11. 

With the assumptions, the governing equations for the liquid 
solution flow were, 


07 Us 

0=n,—> 48 
Hs ay? +ps& (48) 

aT; PT: 
sar (49) 

SOx ts aya 

Xw 8? Xw 

Us Ox Duns ay? (50) 


dP Ug 
0 Taz aye (51) 
uale = De got (52) 
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u X% D PX, (53) 
“az ay3 
At the interface, the mass and energy balances also existed, 
OXw oXy 
psPms oy as —PaDmagy (54) 
dTs, dTa aX, 
— kay, = Kagy, + Padma 4 (55) 
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Fig. 11. Schematic diagram of control volume for a three-dimensional model [26]. 


The above equations can be discretized along the three differ- 
ent directions. Being combined with the boundary and interfacial 
conditions, the equations could be solved. 

Ali et al. [67-69] had used the same mathematical model to 
study the effect of the flow configuration, the inclination angle, the 
Reynolds numbers, various inlet parameters, cu-ultrafine particles 
volume fraction, and thermal dispersion on the performance of the 
dehumidifier/regenerator. 

Mesquita et al. [70] compared three different numerical models 
for parallel plate internally cooled liquid desiccant dehumidifier. 
The second one introduced a constant film thickness. It was 
assumed that the wall was isothermal so that the water flow 
could be neglected. On the basis of some other assumptions, the 
dehumidification process could be described by a two dimensional 
model, as shown in Fig. 12. Being different from the first model 
which does not consider the film thickness, this model took the 
momentum equations into consideration. Firstly, the velocity 
profiles of the air and liquid desiccant were obtained with the 
momentum equations, presented as follows, 


= < (56) 
psS 

Ue = we 50" a(t y)| (58) 

Pa 3p ia | eo 69 


The energy and species equations of the liquid phase, gas 
phase, and the energy and species balances equations at the 
interface are almost the same with those presented in literature 
[26]. Then with the velocity values and boundary conditions, all 
the above equations could be solved numerically and simulta- 
neously with the software package Microsoft EXCEL. 

Recently, Dai and Zhang [71] employed the uniform liquid film 
thickness to evaluate the performance of a cross flow liquid 
desiccant air dehumidifier packed with honeycomb papers. The 
objective of the paper is to analyze the Nusselt and Sherwood 
numbers in the channels with honeycomb papers as the packing 
materials. 


4.3. Models considering variable liquid film thickness 


The final model introduces a variable film thickness. Except the 
constant thickness model, Mesquita et al. [70] also established 


y Symmetry line 


Fig. 12. Schematic diagram of control volume for a two-dimensional model [70]. 
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Table 3 
Detail information of the mathematical models for internally cooled dehumidifier. 


Classification Model 


Flow pattern (air/desiccant) 


Flow pattern (desiccant/cooling fluid) Cooling fluid 


Regardless of film thickness Khan et al. [62] Counter Counter Water and air 
Saman et al. [63] Counter Cross Water and air 
Liu et al. [21] Six different configurations Water 
Yin et al. [25] Cocurrent Cocurrent Water 
Khan et al. [28] Cross Cross Ammonia 
Ren et al. [61] Four possible flow arrangements - 

Uniform film thickness Park et al. [26] Cross Cross R22 
Ali et al. [67-69] Cocurrent/counter/cross - - 
Mesquita et al. [70] Counter Counter Water 
Dai et al. [71] Cross - - 

Variable film thickness Mesquita et al. [70] Counter Counter Water 
Hueffed et al. [72] Cross - - 
Peng et al. [73] Counter - - 
Diaz [74] Cocurrent - - 


another variable thickness model for internally cooled liquid- 
desiccant dehumidifiers. In the model, the film thickness in 
Eq. (56) varied, thus for every step of calculation, the film 
thickness was recalculated, so was the liquid velocity profile in 
Eq. (57). However, to reduce the computational time, the change in 
the airflow velocity profile was neglected as a result of the small 
film thickness changes. Compared the results of variable thickness 
model with the constant thickness model, it was found that the 
two models converged to the same results at higher desiccant flow 
rates. However, the constant thickness model underestimated the 
dehumidification for low desiccant flow rates. 

Hueffed et al. [72] presented a simplified model for a parallel- 
plate dehumidifier, with both adiabatic and isothermal absorption. 
The model used a control volume approach and accounted for 
the film thickness variation by imposing its effect on the heat 
and mass transfer coefficients. The specific equations of the heat 
transfer coefficient in terms of the film thickness were listed as 


Nul 
ac = o (60) 
dp = 2(W — 28) (61) 


Then the mass transfer coefficient was obtained from the 
Chilton-Coulburn analogy as 


-2/3 
ac Di 
= 62 
ap Cp.a (a) (62) 


In each control volume, various parameters, including the film 
thickness ô, hydraulic diameter dą, heat transfer coefficient ac, and 
mass transfer coefficient «p were calculated on the basis of the 
inlet conditions. In this way, the impact of the film thickness 
variation on the heat and mass transfer process were introduced 
into the model. 

Recently, Peng and Pan [73] investigated the transient perfor- 
mance of the liquid desiccant dehumidifier with a one- 
dimensional non-equilibrium heat and mass transfer model. 
Unlike the previous study, the local volume average equations of 
heat and mass transfer were developed in the work, which were 
solved by TriDiagonal-Matrix Algorithm (TDMA). 

In the later year, Diaz [74] also developed a transient two 
dimensional model for a parallel-flow liquid-desiccant dehumidi- 
fier. The geometrical model was almost the same with that of 
Mesquita et al. [70]. The difference of the governing equations lies 
in that the present model took the time item into consideration. 
Some non-dimensional parameters were used to simplify the 
calculation process. With the model, the variations of some critical 


variables over time were illustrated and the effects of oscillatory 
behavior were analyzed in depth. 


4.4. Summary 


The detail information of the mathematical models for the 
internally cooled dehumidifier is listed out in Table 3. 

The models for internally cooled dehumidifier without con- 
sidering the film thickness are very similar to those used in the 
adiabatic dehumidifier. However, the former ones are more 
complicated due to the involvement of the cooling fluid. These 
models ignore the effect of the velocity field. Thus, the results of 
these models have certain discrepancy with the reality, as velocity, 
mass and energy have strong coupling relationship. 

In the models considering uniform film thickness, the film 
thickness and velocity are usually calculated at the beginning 
and then keep constant through the whole calculation. The 
simulation results justified that the constant thickness models 
under-predicted the dehumidification, especially for low desiccant 
mass flow rate. 

For the models considering variable film thickness, all of the 
velocity, mass and energy equations are solved together. In the 
process of each iteration, the film thickness and velocity change, so 
the influence of the flow on the heat and mass transfer can be 
demonstrated. Thus, the model is most accurate. 


5. Conclusion 


Various mathematical models have been proposed to assess the 
performance of the liquid desiccant dehumidifiers. For the adia- 
batic dehumidifier, there are mainly three kinds of models: finite 
difference model, e—NTU model, and simplified models. The finite 
difference model is used most widely for its accuracy. However, it 
involves complicated iterative process which increases the com- 
puter time, so it is only suitable for the design of the component 
and optimization of the operating conditions. Like the finite 
difference model, the e~—NTU model also requires iteration, but it 
is more effective and less accurate correspondingly. The simplified 
models require more assumptions, so they are only suitable in 
certain operation conditions. But due to their high efficiency, they 
are often used to predict the annual performance of the system. 

For the internally cooled dehumidifier, there are also three 
kinds of models: models without considering liquid film thickness, 
models considering uniform liquid film thickness, and models 
considering variable liquid film thickness. The first model ignores 
the effect of the velocity field totally. Thus, the calculation results 
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have certain discrepancy with the reality. The second one consider 
the effect of the velocity field, but the assumption of the steady 
state and fully developed field cannot describe the real condition, 
especially for low desiccant mass flow rate. The final model is most 
accurate as it comprehensively regards the influence of the 
changed velocity field, but the calculation becomes relatively 
complicated. 

Though a large amount of researches have been conducted on 
modeling and analyzing the liquid desiccant dehumidifier, further 
efforts are still needed: 


(1) The vast majority of the previous models assume the heat and 
mass transfer process to be steady state and more transient 
models to study the dynamic performance are needed. 

(2) More three-dimensional models close to the reality are 
needed. 

(3) Most of the models focus on the outlet parameters and more 
research are required to study the heat and mass transfer 
process in the dehumidifier interior. 

(4) It needs to take into consideration the effect of variable 
physical properties which are taken as constant in almost all 
existing models. 
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